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A new method for combining QCD matrix elements and parton showers in Monte Carlo 
simulations of hadronic final states is outlined. The aim is to provide at least a leading-order 
I 1 description of all hard multi-jet configurations together with jet fragmentation to next-to- 

leading logarithmic accuracy, while avoiding the most serious problems of double counting. 



1 Introduction 



The Monte Carlo simulation of multi-jet final states is a challenging problem in QCD and im- 
portant for new physics searches. Two extreme approaches to this problem can be formulated 
as follows. One can use the corresponding matrix elements with bare partons representing jets. 
Then one must add a model for conversion of the partons into hadrons; any realistic model 
will include parton showering, and hence extra jet production and potential double counting. 
Alternatively, one can use the parton model to generate the simplest relevant final state (e.g. 
e + e _ — ► qq) and produce addition jets by parton showering. However, this gives a poor simula- 
tion of configurations with several widely separated jets. 

For earlier work on combining these approaches see refs.Sfl'i&i. Here I outline a method! 
in which the matrix element and parton shower domains are separated at some value y\ of the 
k,T (Durham) jet resolution 13 

y tj = 2mm{ElE 2 j }{\ - cos%)/Q 2 . 

The proposed method has the following features: At > y\ multi-jet cross sections and 
distributions are given by matrix elements modified by Sudakov form factors. At y^ < y\ they 
are given by parton showers subjected to a 'veto' procedure, which cancels the y\ dependence 
of the modified matrix elements to next-to-leading logarithmic (NLL) accuracy. 
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Note that the procedure does not aim at a complete description of any configuration to 
next-to-leading order (NLO) in a s , although this might be possible after subtracting NLO terms 
from the Sudakov form factors (see ref.B). The main objective is to describe all hard multi-jet 
configurations to leading order, i.e. 0{a™~ 2 ) for n jets, together with jet fragmentation to NLL 
accuracy, while avoiding major problems of double counting and/or missed phase-space regions. 



2 Modified Matrix Elements 



The exclusive e + e n-jet fractions at cm. energy Q and fey-resolution Vi = Qi/Q 2 are given to 
NLL accuracy byu 

R 2 (Qi,Q) = [A q (Qi,Q)} 2 

R 3 (Qi,Q) = 2[A,(Q 1 ,Q)] 2 [ Q dqT q (q,Q)A g (Q 1 ,q) 

JQi 

Ra(Qi,Q) = 2[A q (Q 1 ,Q)} 2 \ [ Q dqT q (q,Q)A g (Q 1 ,q) f dq' T q (q', Q)A g (Q 1 , q') 

[JQi JQi 

[ Q dqT q (q,Q)A g (Q 1 ,q) f dq' T g (q\ q)A g (Q 1 , q') 
JQi JQi 



+ 



+ 



Q dqr q (q,Q)A g (Q 1 ,q) f dq' T f {q')A f {Q u q') 
Qi JQi 



etc., where T q ^ g j are q — ► qg, g — > gg and g — > qq branching probabilities 



r. (l ,Q) = (,J _| 

7T q \ q 12 
6TT q 

and A q)g are the quark and gluon Sudakov form factors 



A 3 (Qx,Q) = exp^- J® dqT q (q,Q)^ 
A g (QuQ) = exp(-J*dq\r g (q,Q)+T f (q)] 



with 

Af{Qi,Q) = [A q (Qx,Q)f /A g (Qx,Q) . 

The Sudakov form factor Ai(Qi,Q) represents the probability for a parton of type i to 
evolve from scale Q to scale Q\ without any branching (resolvable at scale Qi). Thus i?2 is the 
probability that the produced quark and antiquark both evolve from Q to Q\ without branching. 
More generally, the probability to evolve from Q to q > Qi without branching (resolvable at 
scale Qi) is A»(Qi, Q)/Ai(Qi, q). 

In i?3, the quark q (or antiquark q) evolves from Q to Q\ without branching, while the q 
(or q) evolves from Q to q, branches, and the resulting partons evolve from q to Q\ without 
branching. 




Thus the overall NLL probability is 

2A q (Q 1 ,Q) Aq f l ' Q) T q (q, Q)A q (Q 1 ,q)A g (Q 1 ,q) , 

which gives Rz{Qi,Q) after integration over Q\ < q < Q. 

We can improve the description of 3-jet distributions throughout the region y qg ,y gg > y\ 
by using the full tree-level matrix element squared -Mqqg in place of the NLL branching prob- 
ability T q (q,Q). We first generate qqg momentum configurations according to M. qqg , with hp- 
resolution cutoff yij > y\ = Q 2 /Q 2 , then weight each configuration with an extra factor of 
[A q (Qi, Q)] 2 A g (Qi, q) where q 2 = mm{y qg , y qg }Q 2 ■ For consistency we also use a s (q) in M. qqg . 

For four or more jets, there are several branching configurations with different colour factors. 
For example there is a contribution from q — > qg branching at scale q followed by g — > gg at 
scale q': 




The probability of this is 

which contributes to the term with colour factor CfCa- The product T q (q,Q)T g (q' , q) is an 
approximation to the full matrix element squared -M qqgg in the kinematic region where y gg is 
smallest interparton separation. Thus it is legitimate in NLLA to replace it by M. qqgg in that 
region. The remaining factor of [A q (Qi, Q)] 2 A g (Qi, q)A g (Qi, q') is the extra Sudakov weight to 
be applied. 

In general, the proposed procedure for generating n-parton configurations is thus as follows: 

• First distribute the parton momenta according to the relevant n-parton matrix element 
squared Ai n , using a fixed value a s (Qi) for the strong coupling. 

• Use the fey-clustering algorithm to determine the resolution values yi = 1 > 2/3 > . . . , > 
y n > yi at which 2,3, ... ,n jets are resolved. These give the nodal values of qj = Q^/yJ 
for a tree diagram that specifies the /cr-clustering sequence for that configuration. 

• Apply a coupling-constant weight factor of a s (qs)a s (qi) • • • a s (q n ) /[a s (Qi)] n ~ 2 < 1. 

• For each internal line of type i from a node at scale qj to qt < qj, apply a Sudakov weight 
factor Ai(Qi, qj) / Ai(Qi, qk) < 1. For an external line from a node at scale qj, the weight 
factor is Ai(Qi,qj). 



Since the weight factors are all less than unity, unweighted events can be generated by rejecting 
those for which the product of weights is less than a random number. 
As an example, the following clustering sequence for e + e~ — > qqggg 




has Sudakov weight 

\(Qi,Q) A q (Qi,Q) A g (Qi, g3 ) 
Ag(<2i, 93) A g (Qi, q 4 ) Ag(Qi, q 5 ) 



x A q (Q 1 ,q 3 )A q (Q 1 ,q 4 )A g (Q 1 ,q 4 ) [A g (Q u q 5 )] 2 
= [A q (Q 1 ,Q)] 2 A g ( y Q 1 ,q 3 )A g {Q u q 4 )A g (Q 1 ,q 5 ) 



Note that the weight factor is actually independent of the structure of the clustering tree and is 
the same as that for the Abelian (QED-like) graph with the same nodal scale values {qj}- The 
clustering of partons will sometimes be 'wrong' but this should not affect LL and NLL terms. 
Other clustering procedures can be envisaged □ which should be equivalent in the dominant 
regions. 



3 Vetoed Parton Showers 



Having generated multijet distributions above the resolution value y\ according to matrix ele- 
ments modified by form factors, it remains to generate distributions at lower values of y^ by 
means of parton showers. This should be done in such a way that the dominant (LL and NLL) 
dependence on the arbitrary parameter y\ cancels. Any residual dependence on y\ could be 
useful for tuning less singular terms to obtain optimal agreement with data. 

Note that y\ must set an upper limit on interparton separations yij generated in the showers. 
Otherwise the exclusive jet rates at resolution y\ could be changed by showering. At first sight, 
this might suggest that we should evolve the showers from the scale Q\ = Qy/yi instead of Q. 
However, this would not lead to cancellation of dependence on logyi. 

Consider, for example, the 2-jet rate at resolution yo = Qo/Q 2 < Ui- If we start from R2 
at scale Q\ and then evolve from Q\ to Qq, we obtain a 2-jet rate of [A q (Qi, Q)A q (Qo, Qi)] 2 
instead of the correct result R2(Qo,Q) = [A q (Qo, Q)] 2 - This is because, although y^ values in 
the showers are limited by y\, the angular regions in which they evolve should still correspond 
to scale Q rather than Q\. Consequently we should allow the showers to evolve from scale Q 
but veto any branching with scale q > Q\ — i.e., the selected parton branching is forbidden but 
that parton has its scale reset to q for subsequent branching. 

The 2-jet rate at any scale Qq < Q\ is now given by the sum of probabilities of 0, 1, 2, . . . 
vetoed branchings (represented by crosses) and no actual resolved branchings: 



The sum of these probabilities for the quark line is 



A q (Q 1 ,Q)A q (Q ,Q)\l+ [ Q dqr g (q,Q) + dqT g (q,Q) f dq' F q (q', Q) + ■ ■ ■ \ 
{ Jqi Jqi Jqi J 

The series sums to 1/ A q (Qi,Q), cancelling the y\ dependence and giving A q (Qo,Q). Similarly 
for the antiquark line. 

For the 3-jet rate at scale Qo < Q\ there are two possibilities: either the event is a 2-jet at 
scale Qi and then has one branching resolved at scale Qo, or it is a 3-jet at scale Qi and remains 
so at scale Qo- The probability of the first case is 




while that of the second case is 



2[A 9 (Q 1 ,Q)] 2 




The sum is indeed 7/1-independent and equal to Rz{Qq, Q)- Similarly for higher jet multiplicities. 

Notice that evolution after a branching at scale q > Qi starts at scale q rather than Q or 
Qi. In general, vetoed showers should evolve in the phase space for angular- ordered branching of 
each partonu. This depends on the colour structure of the matrix element. As illustrated below, 
the angular region for parton i is a cone bounded by the direction of parton j (and vice- versa), 
where i and j are colour-connected. 




If the colour structure is not unique, colour connections must be selected according to their 
relative contributions to the matrix element squared, which are well-defined in the limit that 



the number of colours N c is large. Corrections to the large-iV c limit are normally of relative 
order so this approximation is adequate to ~ 10%. For high parton multiplicity, when the 

colour structure is not known even at large N c , it may be possible to use the clustering scheme 
discussed above as a first approximation in assigning colour connections. 



4 Comments/Conclusions 

• Modified matrix elements plus vetoed parton showers, interfaced at some value y\ of the 
fey-resolution parameter, should provide a convenient way to describe simultaneously the 
hard multi-jet and jet fragmentation regions. 

• The matrix element modifications are coupling-constant and Sudakov weights computed 
directly from the /cT-clustering sequence. 

• Dependence on y\ is cancelled to NLLA by vetoing yij > y\ in the parton showers. 

• This prescription avoids double-counting problems and missed phase-space regions. 

• In principle one needs the matrix elements Ai n for > y\ at all values of n. In practice, 
if we have n < N, then y\ must be chosen large enough for R n >N(Qi, Q) to be negligible. 

• This approach is being implemented (with N = 5) in the e + e~ event generator APACIC++IH 

• It should be possible to extend it to lepton-hadron and hadron-hadron collisions. 

• Extension to NLO along the lines of ref.H may also be possible. 
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